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Abstract 

The helicity density matrix elements poo of p(770) ± and c<j(782) mesons produced in Z° 
decays have been measured using the OPAL detector at LEP. Over the measured meson energy 
range, the values are compatible with 1/3, corresponding to a statistical mix of helicity —1,0 
and +1 states. For the highest accessible scaled energy range 0.3 < xe < 0.6, the measured 
Poo values of the p ± and the u are 0.373 ± 0.052 and 0.142 ± 0.114, respectively. These results 
are compared to measurements of other vector mesons. 



(To be submitted to European Physical Journal C) 



The OPAL Collaboration 



G. Abbiendi 2 , K. Acker staff 8 , G.Alexander 23 , J.Allison 16 , N. Altekamp 5 , K.J.Anderson 9 , 
S.Anderson 12 , S.Arcelli 17 , S.Asai 24 , S.F.Ashby 1 , D. Axen 29 , G. Azuelos 18 < a , A.H.Ball 8 , 

E.Barberio 8 , R.J. Barlow 16 , J.R.Batley 5 , S.Baumann 3 , J. Bechtluft 14 , T.Behnke 27 , 

K.W.Bell 20 , G.Bella 23 , A. Bellerive 9 , S. Bentvelsen 8 , S.Bethke 14 , S.Betts 15 , O.Biebel 14 , 

A.Biguzzi 5 , I.J. Bloodworth 1 , P. Bock 11 , J. Bohme 14 , D. Bonacorsi 2 , M. Boutemeur 33 , 

S.Braibant 8 , P. Bright-Thomas 1 , L. Brigliadori 2 , R.M.Brown 20 , H.J. Burckhart 8 , 

PCapiluppi 2 , R.K.Carnegie 6 , A.A.Carter 13 , J.R.Carter 5 , C.Y.Chang 17 , D.G.Charlton 1 ' 6 , 

D. Chrisman 4 , C. Ciocca 2 , P.E.L. Clarke 15 , E.Clay 15 , I.Cohen 23 , J.E. Conboy 15 , O.C.Cooke 8 , 
J. Couchman 15 , C. Couyoumtzelis 13 , R.L.Coxe 9 , M. Cuffiani 2 , S.Dado 22 , G.M. Dallavalle 2 , 
R.Davis 30 , S.De Jong 12 , A. de Roeck 8 , PDervan 15 , K.Desch 27 , B. Dienes 32 - h , M.S. Dixit 7 , 
J.Dubbert 33 , E. Duchovni 26 , G.Duckeck 33 , LP Duerdoth 16 , PC Estabrooks 6 , E.Etzion 23 , 

F.Fabbri 2 , A.Fanfani 2 , M.Fanti 2 , A.A.Faust 30 , L.Feld 10 , F.Fiedler 27 , M.Fierro 2 , I. Fleck 10 , 

A.Frey 8 , A.Furtjes 8 , D.I. Futyan 16 , P. Gagnon 7 , J.W.Gary 4 , J.Gascon 18 , G. Gaycken 27 , 

C. Geich-Gimbel 3 , G. Giacomelli 2 , P. Giacomelli 2 , V. Gibson 5 , W.R. Gibson 13 , 

D.M. Gingrich 30 ' , D. Glenzinski 9 , J.Goldberg 22 , W. Gorn 4 , C.Grandi 2 , K.Graham 28 , 

E. Gross 26 , J.Grunhaus 23 , M. Gruwe 27 , C.Hajdu 31 G.G.Hanson 12 , M.Hansroul 8 , M.Hapke 13 , 

K. Harder 27 , A. Harel 22 , C.K.Hargrove 7 , M. Harin-Dirac 4 , M. Hauschild 8 , CM. Hawkes 1 , 

R.Hawkings 27 , R.J.Hemingway 6 , CHerten 10 , R.D. Heuer 27 , M.D. Hildreth 8 , J.C.Hill 5 , 

P.R. Hobson 25 , A.Hocker 9 , K.Hoffman 8 , R.J. Homer 1 , A.K. Honma 28 ' a , D. Horvath 31 ' c , 

K.R. Hossain 30 , R.Howard 29 , P. Hiintemeyer 27 , P. Igo-Kemenes 11 , D.C Imrie 25 , K.Ishii 24 , 

F.R.Jacob 20 , A. Jawahery 17 , H.Jeremie 18 , M. Jimack 1 , C.R.Jones 5 , P. Jovanovic 1 , T.R. Junk 6 , 

N. Kanaya 24 , J. Kanzaki 24 , D. Karlen 6 , V. Kartvelishvili 16 , K.Kawagoe 24 , T.Kawamoto 24 , 

PI.Kayal 30 , R.K. Keeler 28 , R.C Kellogg 17 , B.W.Kennedy 20 , D.H.Kim 19 , A.Klier 26 , 

T.Kobayashi 24 , M. Kobel 3 ' d , T.PKokott 3 , M.Kolrep 10 , S. Komamiya 24 , R.V. Kowalewski 28 , 

T.Kress 4 , PKrieger 6 , J. von Krogh 11 , T.Kuhl 3 , PKyberd 13 , CD. Lafferty 16 , H. Landsman 22 , 

D.Lanske 14 , J.Lauber 15 , I. Lawson 28 , J.C Layter 4 , D. Lellouch 26 , J.Letts 12 , L. Levinson 26 , 

R.Liebisch 11 , B.List 8 , C Littlewood 5 , A.W.Lloyd 1 , S.L.Lloyd 13 , F.K. Loebinger 16 , 

G.D.Long 28 , M.J.Losty 7 , J. Lu 29 , J.Ludwig 10 , D.Liu 12 , A. Macchiolo 18 , A. Macpherson 30 , 

W.Mader 3 , M.Mannelli 8 , S. Marcellini 2 , A.J.Martin 13 , J.P.Martin 18 , G.Martinez 17 , 

T. Mashimo 24 , P. Mattig 26 , W.J. McDonald 30 , J. McKenna 29 , E.A. Mckigney 15 , 

T.J. McMahon 1 , R.A. McPherson 28 , F. Meijers 8 , P. Mendez-Lorenzo 33 , F.S. Merritt 9 , H.Mes 7 , 

A. Michelini 2 , S.Mihara 24 , C Mikenberg 26 , D.J.Miller 15 , W.Mohr 10 , A. Montanari 2 , T.Mori 24 , 

K.Nagai 8 , I. Nakamura 24 , H.A. Neal 12 - 3 , R.Nisius 8 , S.W. O'Neale 1 , F.G.Oakham 7 , F. Odorici 2 , 

H.O.Ogren 12 , A.Okpara 11 , M.J. Oreglia 9 , S. Orito 24 , CPasztor 31 , J.R. Pater 16 , G.N.Patrick 20 , 

J.Patt 10 , R.Perez-Ochoa 8 , S.Petzold 27 , P. Pfeifenschneider 14 , J.E. Pilcher 9 , J. Pinfold 30 , 

D.E. Plane 8 , P. Poffenberger 28 , B.Poli 2 , J.Polok 8 , M. Przybycieh 8 - 6 , A. Quadt 8 , CRembser 8 , 

H.Rick 8 , S.Robertson 28 , S.A.Robins 22 , N.Rodning 30 , J.M.Roney 28 , S. Rosati 3 , K.Roscoe 16 , 

A.M.Rossi 2 , Y.Rozen 22 , K.Runge 10 , O. Runolfsson 8 , D.R.Rust 12 , K.Sachs 10 , T.Saeki 24 , 

O.Sahr 33 , W.M.Sang 25 , E.K.G. Sarkisyan 23 , C Sbarra 29 , A.D. Schaile 33 , O.Schaile 33 , 

P. Scharff-Hansen 8 , J. Schieck 11 , S.Schmitt 11 , A. Schoning 8 , M.Schroder 8 , M.Schumacher 3 , 

C. Schwick 8 , W.G. Scott 20 , R. Seuster 14 , T.G. Shears 8 , B.C. Shen 4 , 

CH.Shepherd-Themistocleous 5 , P.Sherwood 15 , G.P. Siroli 2 , A. Sittler 27 , A. Skuja 17 , 

A.M.Smith 8 , G. A. Snow 17 , R.Sobie 28 , S. Soldner-Rembold 10 ^, S. Spagnolo 20 , M. Sproston 20 , 

A. Stahl 3 , K. Stephens 16 , J. Steuerer 27 , K. Stoll 10 , D. Strom 19 , R. Strohmer 33 , B. Surrow 8 , 

S.D.Talbot 1 , P. Taras 18 , S.Tarem 22 , R. Teuscher 9 , M. Thiergen 10 , J.Thomas 15 , 



M.A.Thomson 8 , E. Torrence 8 , S. Towers 6 , I. Trigger 18 , Z. Trocsanyi 32 , E. Tsur 23 , 

M.F. Turner- Watson 1 , I. Ueda 24 , R. Van Kooten 12 , P. Vannerem 10 , M. Verzocchi 8 , H. Voss 3 , 

F.Wackerle 10 , A.Wagner 27 , C.P.Ward 5 , D.R.Ward 5 , P.M. Watkins 1 , A.T.Watson 1 , 

N.K.Watson 1 , P.S.Wells 8 , N.Wermes 3 , D. Wetterling 11 J.S.White 6 , G.W.Wilson 16 , 

J.A.Wilson 1 , T.R.Wyatt 16 , S. Yamashita 24 , V. Zacek 18 , D. Zer-Zion 8 



1 School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 

2 Dipartimento di Fisica dell' Universita di Bologna and INFN, 1-40126 Bologna, Italy 

3 Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 

4 Department of Physics, University of California, Riverside CA 92521, USA 

5 Cavendish Laboratory, Cambridge CB3 OHE, UK 

6 Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ottawa, 

Ontario K1S 5B6, Canada 

7 Centre for Research in Particle Physics, Carleton University, Ottawa, Ontario K1S 5B6, 

Canada 

8 CERN, European Organisation for Particle Physics, CH-1211 Geneva 23, Switzerland 

9 Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 60637, 

USA 

10 Fakultat fur Physik, Albert Ludwigs Universitat, D- 79104 Freiburg, Germany 

11 Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 

12 Indiana University, Department of Physics, Swain Hall West 117, Bloomington IN 47405, 

USA 

13 Queen Mary and Westfield College, University of London, London El 4NS, UK 

14 Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, D-52056 

Aachen, Germany 

15 University College London, London WC1E 6BT, UK 

16 Department of Physics, Schuster Laboratory, The University, Manchester M13 9PL, UK 

17 Department of Physics, University of Maryland, College Park, MD 20742, USA 

18 Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 

Canada 

19 University of Oregon, Department of Physics, Eugene OR 97403, USA 

20 CLRC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK 

22 Department of Physics, Technion- Israel Institute of Technology, Haifa 32000, Israel 

23 Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 

24 International Centre for Elementary Particle Physics and Department of Physics, University 

of Tokyo, Tokyo 113-0033, and Kobe University, Kobe 657-8501, Japan 

25 Institute of Physical and Environmental Sciences, Brunei University, Uxbridge, Middlesex 

UB8 3PH, UK 

26 Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 

27 Universitat Hamburg/DESY, II Institut fur Experimental Physik, Notkestrasse 85, D-22607 

Hamburg, Germany 

28 University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, Canada 

29 University of British Columbia, Department of Physics, Vancouver BC V6T 1Z1, Canada 

30 University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 

31 Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, Hungary 

32 Institute of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 



33 Ludwigs-Maximilians-Universitat Miinchen, Sektion Physik, Am Coulombwall 1, D-85748 
G arching, Germany 



a and at TRIUMF, Vancouver, Canada V6T 2A3 

6 and Royal Society University Research Fellow 

c and Institute of Nuclear Research, Debrecen, Hungary 

d on leave of absence from the University of Freiburg 

e and University of Mining and Metallurgy, Cracow 

' and Heisenberg Fellow 

9 now at Yale University, Dept of Physics, New Haven, USA 

h and Depart of Experimental Physics, Lajos Kossuth University, Debrecen, Hungary. 



1 Introduction 



Very little is known about the role of spin in the hadronization process. At LEP, this can 
be investigated by studying the properties of vector mesons produced in hadronic Z° decays. 
Recent data on the helicity states of vector mesons produced in hadronic Z° decays |l|-§| reveal 
that the spin of high-energy K*(892)° and 0(1020) mesons is preferentially aligned transverse 
to the direction of their momentum. Measurements of the p(770)° and D*(2010) ± mesons 
are consistent with this behaviour, while B* mesons show no alignment. The mechanism at 
the origin of spin alignment is not well understood theoretically |J, and this phenomenon is 
ignored in models such as the Lund string model || and the cluster model [0]. Extending 
these studies to other vector mesons would provide valuable help in elucidating the role of spin 
in hadronization. The helicity density matrix elements for the p(770) ± and o>(782) mesons 
produced in Z° decays have never been measured: these are particularly interesting because 
the most pronounced alignments observed so far are for light mesons. The p 1 * 1 and u mesons, 
together with the p° whose helicity matrix element has already been measured 0, have a 
similar quark structure and can be expected to have similar behaviour. Experimentally, the 
systematic uncertainties affecting the extraction of the helicity matrix elements for these three 
mesons differ substantially. 

This paper describes the measurement of the helicity density matrix elements poo of p ± 
and us mesons produced in Z° decays using the OPAL detector at LEP. The detector and its 
performance are described in detail in Refs. |8HlQf. The data sample consists of 4.1 million 



hadronic Z decays collected at centre-of-mass energies within ±2 GeV of the Z peak. The 



selection of hadronic events is presented in Ref . |Tl[ . The method to reconstruct and identify 
the p ± — > ir^iT and u — > tt+tt^tt decays is the same as used in Ref. |)12|| , where it is explained 
in detail. The measurements of poo for the p^ and u mesons are presented in Sections 2 and 3, 
and the results are compared to those for other mesons at LEP in Section 4. 



2 Helicity density matrix element poo °f P mesons 

The p ± resonance decays dominantly via the ^tt channel. Using as a spin analyser the angle 
in the tt^ti rest frame between one of the pion momenta and the p^ boost direction, the 
distribution of this angle, 9h, is [jl3| : 



W(cos6 H ) = |j[(l-poo) + (3poo-l)cos 2 ^] (1) 

where poo is the helicity density matrix element expressing the probability that the spin of the 
p ± meson be perpendicular to its momentum direction. 

As the mechanisms producing the observed mesons and their amount of spin alignment 
can vary with their energy, the analysis is repeated for different intervals of the scaled energy 

•&E -C-'meson / -'-'beam • 



2.1 p^ meson reconstruction 

The reconstruction of the decay p ± — > 7r ± 7r° follows exactly that described in Ref. [12]. Charged 
pion candidates are selected as tracks in the central drift chambers with an energy loss mea- 



surement having a probability greater than 1% for the pion hypothesis | 1CJ . Neutral pions are 
obtained from the combination of pairs of photons detected either as localised energy deposits 
in the electromagnetic calorimeter or as two tracks from a 7 — > e + e~ conversion within the 
volume of the central drift chambers, and selected using the multivariate method described in 
Ref. 0. 



All selected 7r° and tc^ candidates are combined in pairs. From the energy and momenta 
of the 7T° and tt^ in each pair, three quantities are evaluated: the scaled energy of the pair, 
xe, the invariant mass of the tt^tt system, m, and the cosine of the spin analyser angle Oh, 
defined as the angle in the 7r ± 7r° rest frame between the 7r° momentum and the boost of the 
n^rr system. The sample is divided into ten equal bins of cos Oh in the range from —1 to 
+ 1, for six intervals of xe in the range from 0.025 to 1. The size of the bins in xe and cosOh 
are large compared to the experimental resolutions on these quantities, which are dominated 
by the ir° energy resolution. The Monte Carlo simulations provide realistic estimates of these 
resolutions as they reproduce well the width of the tt° mass peak ||12|| . They predict that the 
7T° energy resolution varies between 4% and 8% over the energy range relevant for the present 
analysis, resulting in a resolution of approximately 0.04 on cosOh- 

The number of reconstructed p ± mesons in each cos Oh bin and xe interval is obtained from 



a fit to the corresponding invariant mass distribution using the method described in Ref. |L2 
These numbers are corrected by the efficiency evaluated by applying the same reconstruction 
method and fit procedure to a sample of 6.4 million hadronic Z° decays simulated |14| using 



the Monte Carlo programs JETSET 7.3 and 7.4 B tuned to reproduce the global features of 



hadronic events as observed by OPAL [15, 16 . 



Examples of invariant mass distributions and fit results are shown in Figs. [I] and Q. As can 
be seen, the shapes of the p ± signal and of the underlying background vary significantly as a 
function of xe and cosOh- According to the Monte Carlo simulations, this is due in large part to 
the dependence of the 7r° reconstruction efficiency on its energy. To track these effects carefully 
the fit procedure, described in Ref. |l!| and summarised below, is applied independently to 



each mass distribution. 

In the fit, the p^ signal is parameterized as a relativistic Breit-Wigner convoluted with the 

2 2 

experimental mass resolution and multiplied by a factor 1 + C m ° ™ where V is the width 
of the resonance |T7| and m is the Breit-Wigner pole mass. This factor has been shown to 



take adequately into account shape distortions due to interferences and residual Bose-Einstein 
correlations [pl|. The systematic variations of the signal shape consist of fixing C to zero or 
leaving it as a free parameter and, in addition, the mass resolution is set to the Monte Carlo 
prediction or left as a free parameter. The systematics associated with the background are 
evaluated using two methods. In the first, the background shape is taken from the simulation 
and normalised to the number of counts outside the signal region. In the second method, the 
background is parameterized as: 

f(m) = p^m - m th ) P2 x exp(p 3 (m - m th ) + p 4 (m - m th ) 2 ) , (2) 

5 



where rrith = ra^± + m^o and the parameters p\ to p^ are determined in the fits to the data. 
A Gaussian representing the reflection from u — > 7r + 7r _ 7r° decays is added to this shape. Its 
width is fixed to the Monte Carlo prediction while its amplitude and centroid are left as free 
parameters in order to absorb possible imperfections in the modelling of the background near 
the 7T ± 7r° threshold. The ir° selection is also varied as in Ref. [O], testing the sensitivity of the 



results to an increase and a decrease of the acceptance by a factor 2. 



2.2 Extraction of the matrix element p 00 of p ± mesons 

For each interval of xe, the efficiency-corrected p yields are evaluated for all ten cos 8h bins 
and fitted to the expression 

I(cos6 H ) = A(l + Bcos 2 9 H ) , (3) 

From which the value of poo is obtained via: 

1 + B m 

p0 ° = sT^ • (4) 

The parameters A and B and their errors are obtained from a linear least-squares fit. As an 
example, Fig. ^| shows the fits to the data of Figs. [I] and 0. The fit of Eqn. [3] to the data is 
repeated for all systematic variations of the signal and background parameterization, of the ir° 
selection, and of the Monte Carlo sample used for the efficiency. The resulting B values are 
averaged and the systematic error associated to each source of uncertainty is taken from the 
rms deviation from the average. In this way, these errors reflect the uncertainty on B and are 
independent of global variations which affect only the parameter A. 

The measured p o values are listed in Table [3], together with the statistical errors and the 
errors from the following systematic uncertainties: 

1. The statistical error on the Monte Carlo samples used to calculate the efficiency. 

2. The bias induced by the presence of the background under the signal peak. This is 
estimated by the difference of the results of fits to the invariant mass spectra of Monte 
Carlo samples where the background is included or excluded. 

3. The variations of the fitted yields when the parameterization of the shape of the signal 



and the background is varied as in Ref. | 12|| . 



4. The difference in efficiency obtained with Monte Carlo samples using the JETSET tune 
parameters of Refs. |D| and ||16|| . 



5. The variation observed when the analysis is repeated with different values for the cut on 
the 7T° selection variable |i~2| , corresponding to changes by factors from 1/2 to 2 in the 
acceptance. As shown in Ref. []12] , the variation of this cut induces significant changes in 



the shape of the background and therefore provides an additional test of the stability of 
the results relative to the assumptions regarding its shape. 



3 Helicity density matrix element poo °f w mesons 



The decay uo — > 7r + 7r~7r° has a branching ratio of 88.8±0.7% []17 |. In the rest frame of the 
7T + 7r _ 7r system, the momenta of the three pions are in a plane. The appropriate spin analyser, 
9jj, in this case is the angle between the normal to this plane and the boost direction ||13| , and 
Eqn. [l] applies here [19| also. 



3.1 lj meson reconstruction 



The reconstruction of the decay uj — > 7t + 7t~tt follows exactly that described in Ref. [12 
Charged and neutral pion candidates are selected as in the p analysis. All triplets comprising 
two oppositely charged pions and one neutral pion are considered. From the energy and mo- 
menta of the 7r + , 7T~ and tt° candidates three quantities are evaluated: the scaled energy of the 
triplet, xe, the invariant mass of the tt + 7t~7t° system, m, and the cosine of the spin analyser 
angle, Oh, defined above. The sign of cos Oh is arbitrary and only its absolute value is con- 
sidered. The Monte Carlo simulations predict that the resolution on | cos6^| is approximately 
0.04, increasing to 0.06 for 0.3 < xe < 0.6. The sample is divided into six equal bins of | cos0h\ 
in the range from to 1, for six intervals of Xe in the range from 0.025 to 1. 

The number of reconstructed uo mesons in each | cos Oh | bin and xe interval is obtained from 
a fit to the corresponding invariant mass distribution. Examples of invariant mass distributions 
and fit results are shown in Figs. |4] and |5|. In these fits, the combinatorial background is 
parameterized as a third order polynomial, P^{m). As in Ref. [12], the uo signal, S(m), is 



described as the superposition of two Gaussians sharing the same centroid, the ratios of the 
two widths and areas being determined from the simulation. According to the simulations, 
the position and total width of the peak does not depend on [cos^l- Therefore, for each 
Xe interval, these two quantities are determined from a fit to the data in the total interval 
< | cos0 H \ < 1, and are subsequently fixed to these values in the fits to the six individual 
| cos Oh I bins. 



3.2 Extraction of the matrix element p 00 of uo mesons 

The extracted uj yields are corrected by the reconstruction efficiency evaluated in a similar 
manner to that in the p analysis (Section |2.2j ). For each interval in xe, the efficiency-corrected 
uj yields are evaluated for the six | cos6^| bins and fitted to Eqn. |3|. As an example, Fig. |] 
shows the fits to the data of Figs. ^ and |5|. The shape of the background does not vary with 
xe as much as in the case of the p , and here a fit to the sum of all xe bins (Fig. |]) can be 
performed. This allows a high-statistics test of the fitting procedure. The measured poo values 
are listed in Table 0, together with the statistical and systematic uncertainties. The systematic 



errors are evaluated as in Section |2]2j. In the case of the uj, the systematic variation of the shape 
of the signal (third error in table 0) are derived from determining the width of the u> peak from 
the simulation instead of the data. As in Ref. | 12| , an additional error (the sixth in table |2]) 



is derived from the test of whether the extracted uj signal has the expected dependence on the 



matrix element of the uo decay [ED 



t* 1 2 



|pl x p;r 

A w — l - ^ ^775 1-3J 

|p* x pl\ 2 



I m«.r 



where ptji is the momentum of the 7r ± meson in the 7r + 7r~7T° rest frame. The A^ distribution is ex- 
pected to rise linearly with X w for the signal, and be flat for the combinatorial background. The 
data are therefore divided into six bins in A w , and the two-dimensional distributions I(m, XJ) 
are fitted with the expression: 

I{m,\u) = -, rTTn~T^ 6 

e(m,Aa,) (l + Awd) 

where e(m, A w ) is the acceptance of the detector determined from the combinatorial background 
distributions in the simulations, and S is a free parameter in the fit, designed to account for 
possible deviations between the A^-dependence of the acceptance in the data and in the Monte 
Carlo. The value of the sixth error in Table § is the difference between the results obtained 
from the fits to the I(m, A w ) and I(m) distributions, where I(m) is the sum of the six A w bins. 



4 Results and discussion 

Tables [l] and || present the measured poo values for the p ± and uo mesons, respectively, together 
with their statistical, systematic and total errors. The total numbers of meson candidates are 
also listed, together with their statistical errors. The statistical errors are dominated by the 
contribution of the background under the p ± and 00 signals. The highest accessible xe interval 
is 0.3 < Xe < 0.6. The uncertainty on the p^ and uo samples with xe > 0.6 (approximately 
400 ± 200 and 800 ± 400 candidates, respectively) is too large for a meaningful extraction of 



Poo- The systematic errors are evaluated using the method of Ref. |L2|. In Tables 1 and 2, 
the scatter of the error values as a function of xe for some uncertainty sources indicate that 
their evaluation may still be affected by statistical fluctuations. However the observed scatter is 
smaller than the size of the total errors and thus should not affect significantly the final results. 

The measured values of poo as a function of xe are shown in Figs, ^a and b for the p ± 
and uo mesons, respectively. The measurements are compatible with 1/3 (the dashed line in 
Figs 0a and b), corresponding to a statistical mix of helicity —1,0 and +1 states. Out of ten 
measurements, eight are equal to 1/3 within one standard deviation, and two are within two 
standard deviations. In other cases where vector meson spin alignment has been observed |]-[t| 
it appears only at meson energies above xe > 0.3. Here, in the energy range 0.3 < xe < 0.6, 
the values of p 00 for the p ± and u mesons are 0.373 ± 0.035 ± 0.038 and 0.142 ± 0.081 ± 0.080. 

According to the Monte Carlo, and depending on the p energy, up to 10% of p 1 * 1 originate 
from the sequential decays of J p = 0~ mesons into a p ± ( J p = 1~) and another J p = 0~ meson. 
The most important source of these decays is D° — > p + K~. In these decays the orbital angular 
momentum of the p must be opposite to its intrinsic spin so that the p 1 * 1 must be in a helicity 
state (poo = 1) in the rest frame of the parent meson. The degree of alignment in the laboratory 
frame depends on the relative momenta of the p 1 * 1 and the parent 0~ meson. Simulations predict 
that only 60% of the alignment survives in the laboratory frame for p ± mesons in the interval 

8 



0.3 < xe < 0.6 and that it essentially disappears for xe < 0.1. The effect of the alignment 
of the p ± mesons coming from J p = 0~ — > 0~ + 1~ decays present in the JETSET 7.3 and 
7.4 samples based on the parameters of Ref. |T5j and [|16j are almost indistinguishable, despite 



the inclusion of L — 1 mesons in the latter simulation. This prediction is shown in Fig. [FJa 
as a dotted line. It agrees with the data. However, these small deviations from 1/3 at high 
Xe are comparable in size to the experimental errors. For the u, the simulations predict that 
less than 5% of uj mesons in the interval 0.3 < xe < 0.6 come from J p = 0~ — > 0~ + 1~ 
decays. The expected impact on the average u alignment, shown as a dotted line in Fig. 0b, is 
significantly smaller than in the p ± case (Fig. [7(a). If the contributions of the expected sources 
of J p = 0~ — > 0~ + 1~ decays are removed, the remaining population is still compatible with 
Poo — 1/3 over the entire xe range. 

Fig. 0c compares the p 1 * 1 and uo measurements with those for the K*(892)° obtained by 
OPAL H and the p° by DELPHI [|. The less precise DELPHI K*(892)° data @ are not 
shown for clarity: they are consistent with the OPAL data. The K*(892)° data show a significant 
preference for poo values above 1/3 at values of xe above 0.3. By itself, the p° data was not 
sufficiently precise to conclude whether similar alignment values were also observed for light, 
unflavoured mesons in the corresponding xe range. 

Fig. H shows a compilation of poo measurements for different mesons in different xe ranges. 
Up to now, the B* meson was the only case where poo=l/3 was clearly preferred. The K*(892)° 
and <p mesonsQ appeared to prefer larger values of poo > 1/2 at high xe, with the D* mesons in 
between. In contrast, the new results on the p ± and u, together with the previous DELPHI p° 
results, appear to prefer p 00 values close to 1/3. Therefore the presence of spin alignment above 
xe > 0.3 cannot be considered to be a general property of mesons produced in hadronic Z° 
decays. This could either be due to the influence of cascade decays on the observed alignments 
or to some unknown mechanism producing the alignment. More measurements above xe = 0.3, 
and particularly above 0.6, would contribute significantly to the understanding of meson spin 
alignment. 



5 Conclusion 

The helicity density matrix elements poo of p^ and u mesons produced in Z° decays have been 
measured using the OPAL detector at LEP. Over the entire energy range, the measured values 
are compatible with 1/3 corresponding to a statistical mix of helicity —1,0 and +1 states. The 
measurements in the highest accessible energy range 0.3 < xe < 0.6 are 0.373 ± 0.052 and 
0.142 ± 0.114 for p 1 * 1 and u mesons, respectively. Taken together, these results are lower than 
the values observed at high xe for the K*(892)° and </> mesons. 



lr The OPAL </> data ||] are corrected for the effects of J p = ^0 +1 decays predicted by the Monte 
Carlo simulations. However the size of the correction is small compared to the observed alignment. 
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Stat. 


Bkg. 


Signal 


Diff. 


7T° 


Total 


range 


(xlO 3 ) 


value 


(data) 


(MC) 


bias 


shape 


MC 


select. 


error 


0.025 - 0.050 


48 ± 5 


0.312 


0.066 


0.004 


0.000 


0.000 


0.005 


0.047 


0.081 


0.050- 0.100 


156 ± 5 


0.338 


0.060 


0.030 


0.030 


0.008 


0.004 
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0.081 


0.100 - 0.150 


222 ± 4 
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0.008 


0.018 


0.058 


0.150 - 0.300 


273 ± 3 
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0.027 


0.019 
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0.018 
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0.300 - 0.600 


60 ± 1 


0.373 


0.035 


0.029 


0.001 


0.011 


0.009 


0.022 


0.052 



Table 1: Measured poo values for p^ mesons as a function of xe, together with the statistical, systematic and 
total errors. The different systematic error contributions 1-5 are described in Sect. 2.2. The total numbers 
of p^ candidates in each xe bin are also listed, together with their statistical errors, for the 7r° cut yielding 
the largest acceptance. These numbers are only shown to give an idea of the statistical precision of the 
measurement and are not used in the analysis. 
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100.2 ± 1.7 


0.249 


0.018 


0.003 


0.047 


0.008 


0.008 


0.038 


0.007 


0.065 


0.100 - 0.150 


47.5 ± 1.2 


0.308 


0.021 


0.003 


0.015 


0.003 


0.010 


0.032 


0.005 


0.043 


0.150 - 0.300 


41.4 ± 1.0 
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Table 2: Measured poo values for uu mesons as a function of xe, together with the statistical, systematic and 
total errors. The different systematic error contributions 1-6 are described in Sect. 2.2 and 3.2. The total 
numbers of uj candidates in each xe bin are also listed, together with their statistical errors, for the tt° cut 
yielding the largest acceptance. These numbers are only shown to give an idea of the statistical precision of 
the measurement and are not used in the analysis. 
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Figure 1: Invariant mass distributions of tt^tt combinations for the range 0.10 < Xe < 0.15, 
for the ten equal bins of cos6h between —1 and +1. The black circles are the data. The full 
lines are the background obtained in the fits to the data. The white circles show the /r* 1 signal, 
scaled by a factor 2, obtained by subtracting the fitted background (full lines) from the data. 
The dashed lines are the fitted signals scaled by the same factor. 
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Figure 2: Invariant mass distributions of 7r ± 7r combinations for the scaled energy range 0.3 < 
xe < 0.6, for the ten equal bins of cos6h between —1 and +1. The black circles are the data. 
The full lines are the background obtained in the fits to the data. The white circles show the 
p^ 1 signal, scaled by a factor 2, obtained by subtracting the fitted background (full lines) from 
the data. The dashed lines are the fitted signals scaled by the same factor. 
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Figure 3: Efficiency-corrected p ± yields obtained from the fits to the data shown (a) in Fig. 1 
(0.10 < x E < 0.15) and (b) in Fig. 2 (0.3 < x E < 0.6). For clarity, the total yields are 
normalised to 1 using the results of the fit to the data in the entire cos 9h range. The errors 
are statistical only. The full lines represent the fit of Eqn. [3] to the data. 
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Figure 4: Invariant mass distributions of 7r + 7r~7r° combinations for the range Xe > 0.025, for 
the six equal bins of | cos 0jj\ between and 1. The black circles are the data. The full lines 
are the background obtained in the fits to the data. The white circles show the u signal, scaled 
by a factor 10, obtained by subtracting the fitted background (full lines) from the data. The 
dashed lines are the fitted signals scaled by the same factor. 



17 



OPAL (0.3<x F <0.6) 



10000 



> 



10000 



o 
a 

© 



10000 




^^^H^iH 



7i 7C + 7u" Invariant mass (GeV/c ) 



Figure 5: Invariant mass distributions of 7r + 7r~7r° combinations for the scaled energy range 0.3 
< xe < 0.6, for the six equal bins of | cos##| between and 1. The black circles are the data. 
The full lines are the background obtained in the fits to the data. The white circles show the 
uj signal, scaled by a factor 10, obtained by subtracting the fitted background (full lines) from 
the data. The dashed lines are the fitted signals scaled by the same factor. 
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Figure 6: Efficiency-corrected u yields obtained from the fits to the data shown (a) in Fig. 5 
(xe > 0.025) and (b) in Fig. 6 (0.3 < xe < 0.6). For clarity, the total yields are normalised to 
1 using the results of the fit to the data in the entire | cos 9h\ range. The errors are statistical 
only. The full lines represent the fit of Eqn. to the data. 
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Figure 7: Measured p 00 values as a function of xe for (a) p ± mesons and (b) uj mesons produced 
in Z° decay. The dashed lines correspond to an isotropic spin distribution (1/3). The dotted 
lines represent the amount of alignment expected from the J p = 0~ — > _ + 1~ decays present 
in the simulations described in Sect. 4. In (c), the data from (a) and (b) are compared with 
other measurements 0,131. 
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Figure 8: Summary of published poo measurements for vector mesons produced in Z° decays. 
The results for the p and uj are from this work. The Delphi results for the p°, K*(892)° and </> 
mesons are from Ref. ||, and the OPAL results for the <j), D* ± and B* mesons are from Ref. |T]|. 
Other results are from Refs. § (OPAL K*(892)°) and [| (ALEPH and DELPHI B*). Only 
measurements at large x or corresponding to primary quarks are shown. The K*(892)° data at 
0.3 < xe < 0.5 and xe > 0.7 are also shown as an example of the possible xe dependence. 
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